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ABSTRACT 

Eleven candidate damped Lya absorption systems were identified in twenty-seven 
spectra of the quasars from the APM z > 4 survey covering the redshift range 
2.8 < Zabsorption < 4.4 (8 with z a b S orption > 3.5). High resolution echelle spectra 
(0.8A FWHM) have been obtained for three quasars, including two of the highest 
redshift objects in the survey. Two damped systems have confirmed HI column den- 
sities of A^hi^ 10 20 3 atoms cm" 2 , with a third falling just below this threshold. We 
have discovered the highest redshift damped Lya absorber known at z=4.383 in QSO 
BR1202-0725. 

The APM QSOs provide a substantial increase in the redshift path available for 
damped surveys for z>3. We combine this high redshift sample with other quasar sam- 
ples covering the redshift range 0.008 < z < 4.7 to study the redshift evolution and the 
column density distribution function for absorbers with log JVhi> 17.2. In the HI col- 
umn density distribution f(N) = kN~@ we find evidence for breaks in the power law, 
flattening for 17.2 < logTVjji H 21 and steepening for logAni> 21.2. The breaks are more 
pronounced at higher redshift. The column density distribution function for the data 
with log A r ni>20.3 is better fit with the form f(N) = (f*/N*)(N/N*)-Pexp(-N/N*) 
with log N, = 21.63±0.35, (3 = 1.48±0.30, and /* = 1.77 x 10" 2 . We have studied the 
evolution of the number density per unit redshift of the damped systems by fitting the 
sample with the customary power law N(z) = N${\ + z) 1 . For a population with no 
intrinsic evolution in the product of the absorption cross-section and comoving spatial 
number density this will give 7 = 1/2 (O = 1) or 7 = 1 (O = 0). The best maximum 
likelihood fit for a single power law is 7 = 1.3 ± 0.5 and No = .04+q 3 ,, consistent with 
no intrinsic evolution even though the value of 7 is also consistent with that found for 
the Lyman limit systems where evolution is detected at a significant level. However, 
redshift evolution is evident in the higher column density systems with an apparent 
decline in N(z) for z>3.5. 

Key words: cosmology — galaxies: evolution — galaxies: formation — quasars: absorp- 
tion lines— quasars: individual (BR1033-0327, BRI1108-0747, BR1202-0725) 



1 INTRODUCTION 

This paper is the third of a series presenting results from 
studies of the QSOs discovered in the APM survey for z 
> 4 quasars. A study of the evolution of Lyman limit ab- 
sorption systems over the redshift range 0.04 < z < 4.7 was 
presented in Storrie-Lombardi et al. (1994) [Paper I]. The 
intermediate resolution (5A) QSO spectra and the survey 



for high redshift damped Lya absorbers are presented in 
Storrie-Lombardi et al. (1996) [Paper II]. The evolution of 
the cosmological mass density of neutral gas at high redshift 
and the implications for galaxy formation theories are dis- 
cussed in Storrie-Lombardi, McMahon & Irwin (1996) [Pa- 
per IV]. In separate papers we will describe the intrinsic 
properties of the QSOs and studies of the Lya forest clouds 
at high redshift. A high resolution study of the Lya forest 
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region in a redshift z = 4.5 QSO has been completed by 
Williger et al. (1994). 

How and when galaxies formed are questions at the fore- 
front of work in observational cosmology. Absorption sys- 
tems detected in quasar spectra provide the means to study 
galaxy formation and evolution up to redshifts of approxi- 
mately five, back to when the Universe was less than 10 per- 
cent of its present age. Surveys for absorption features have 
several advantages over trying to directly detect galaxies at 
high redshift. Much shorter exposure times are required be- 
cause the QSOs are relatively bright (R~ 18-19.5) and the 
large equivalent width systems are easily detected in the 
spectra. This provides good absorption candidates to follow 
up with higher resolution spectra. The redshift and column 
density can be accurately determined from the wavelength 
of the absorption system and the line profile. This is far eas- 
ier and more reliable than trying to directly get a spectrum 
of a very faint high redshift galaxy. 

While the baryonic content of spiral galaxies that are 
observed in the present epoch is concentrated in stars, in 
the past this must have been in the form of gas. The princi- 
pal gaseous component in spirals is neutral hydrogen which 
has led to surveys for absorbers detected by the damped 
Lya lines they produce (Wolfe et al. 1986, hereafter WTSC; 
Lanzetta et al. 1991, hereafter LWTLMH; Lanzetta, Wolfe & 
Turnshek 1995, hereafter LWT; Wolfe et al. 1995; Paper II). 
Though damped Lya systems are observationally very rare 
objects with ~40 confirmed examples known, the HI mass 
per unit comoving volume they contain is roughly compa- 
rable to the mass density of baryonic matter in present-day 
spirals, ie. a major constituent of the Universe (Wolfe 1987, 
LWT). Their metal abundances are much lower than Galac- 
tic values (Pettini, Boksenberg & Hunstead 1990; Rauch et 
al. 1990; Pettini et al. 1994) and they are characterised by 
low molecular content and low, but not negligible, dust con- 
tent (Fall, Pei & McMahon 1989; Pei, Fall & Bechtold 1991; 
Pettini et al. 1994), features consistent with an early phase 
of galactic evolution. They may be the progenitors of spi- 
ral galaxies like our own and are clearly important for the 
study of the formation and evolution of galaxies. They have 
been detected across a very large redshift range z~ [0.5,4.5] 
providing the means to pinpoint the epoch of formation of 
disk galaxies and study their evolution. 

Eleven candidate damped Lya absorption systems out 
of 32 measured Lya features were identified in 27 spectra of 
the mainly non-BAL quasars from the APM z > 4 survey 
(Paper II). The eleven candidates cover the redshift range 
2.8 < Zatsorption < 4.4 (8 with z a bsor P tion > 3.5) and have 
estimated column densities JVhi> 10 20 ' 3 atoms cm -2 . In this 
paper the QSO BR1 144— 0723 with a candidate absorber at 
z=3.26 is removed from further consideration in the sam- 
ple. It has been observed with the Anglo- Australian Tele- 
scope at high resolution and the damped candidate has been 
found to be all OVI absorption at z=4.0 (R. Hunstead, pri- 
vate communication). High resolution echelle spectra (0.8A 
FWHM) were obtained by S. D'Odorico as part of the ESO 
key programme studying high redshift quasars for four of the 
QSOs in the APM sample (BRI0952-0115, BR1033-0327, 
BRI1108-0747, BR1202-0725). The signal-to-noise ratio in 
BRI0952— 0115 was very poor but the other spectra have 
been used to confirm two Lya features as damped with an- 
other falling just below the log JVhi> 20.3 threshold. We 



have discovered the highest redshift damped Lya absorber 
known at z=4.383 in QSO BR1202-0725. The confirmation 
of the absorption systems is discussed in section 2. These 
data have been combined with data from previous surveys 
(WTSC, LWTLMH, and LWT) and the results for the Ly- 
man limit systems obtained in Paper I to study the HI col- 
umn density distribution for logiVHi> 17.2 and redshift evo- 
lution of these systems for 0.008 < z < 4.7. 

Numerous authors have studied the distribution of col- 
umn densities, f{Nm), for Lya absorption lines. The first 
determination was by Carswell et al. (1984) for lines with 
10 13 < Ahi< 10 16 atoms cm" 2 . They found f(N) cx N~ 
(P — 1.7 ± 0.1). Damped Lya absorption (DLA) systems 
comprise the high column density tail of neutral hydrogen 
absorbers with column densities of Ahi> 2 x 10 20 atoms 
cm~ 2 . They dominate the baryonic mass contributed by HI. 
When damped systems are included in the column density 
distribution function for a single power law fit the expo- 
nent is P = 1.4-1.7 (cf. Tytler 1987; Petitjean et al. 1993 
and references therein) . Assuming the baryonic mass is pro- 
portional to the HI column density and takes the form 
f{Nm) oc N' 13 for the HI column density distribution func- 
tion, the mass contribution from the damped systems can 
be estimated as 

Mtotal oc / Nmf(N m )dNm 

J Ni 

r-N 2 

1 



r«2 

/ N'^NdN (assume (3 ^ 2) 

J JVi 



2-/3 



N 2-0 _ M 2-f3 



(1) 



One problem with the power law representation is that if 
P < 2, as all current estimates indicate, then the total mass 
in damped systems diverges unless an upper bound to the 
HI column density is assumed. For example, if we take 20.3 
S: log Nm ~ 22, the fractional contribution to the total 
HI mass for damped systems, Mf, is then Mf — 0.86 for 
P — 1.5 and Mf — 0.69 for P = 1.7. However, there is no a 
priori reason for assuming this upper limit and hence there 
is no strict upper bound to any estimate of the total HI 
mass in damped systems. An alternative parameterisation 
using a gamma function to describe the HI c olumn density 
distribution was adopted by Pei & Fall (1995) and provides 
an elegant solution to the diverging mass problem. We dis- 
cuss these points in more detail in section 3 and the redshift 
evolution of the absorbers in section 4. 



2 CONFIRMATION OF DAMPED Lya 
SYSTEMS 

2.1 Echelle Observations 

Echelle spectra of four QSOs were obtained in March, 1993 
by S. D'Odorico as part of an ESO key programme studying 
high redshift quasars. They were taken at La Silla with the 
3.5m NTT telescope using the EMMI instrument in echelle 
mode using a 2048x2048 pixel LORAL CCD as the detector. 
A slit of 15" in length was used and generally the slit width 
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Table 1. ESO Observations March, 1993 



Table 2. Lyman-a Absorption Systems 



QSO 




Date (UT) 


Exp 
(sees) 


Grating 


Slit 
(") 


QSO 




Absorption 
Redshift 


log jV HI 


BRI0952 


— 0115 


93 Mar 14 


5400 


GR9CD3 


1.2 


BR1033 — 


0327 


4.14945 


19.80 ± 0.10 






93 Mar 14 


7200 


GR9CD3 


2.0 






4.15314 


18.69 ± 0.30 


BR 1033 


- 0327 


93 Mar 15 


7200, 8000 


GR9CD3 


1.2 






4.16647 


19.70 ± 0.15 


BRI1108 


- 0747 


93 Mar 15 


6000 


GR9CD3 


1.2 






4.16726 


19.37 ± 0.42 






93 Mar 15 


6000 


GR9CD4 


1.2 






4.17481 


19.60 ± 0.24 


BR 1202 


- 0725 


93 Mar 14 


8000 


GR9CD3 


1.2 














93 Mar 15 


6797 


GR9CD3 


1.2 


BRI1108 - 


- 0747 


3.60673 


20.33 ± 0.15 






93 Mar 15 


7200, 8000 


GR9CD4 


1.2 






















BR1202 - 


0725 


4.38290 


20.49 ± 0.15 















was 1.2". Two grating setups were used, one covering 4700- 
8300A and the other covering 5800-9500A with a resolution 
of ~40 km s" 1 (lA). [See Giallongo et al. 1994 for more 
details.] The observations are summarised in table |l[ 

The spectra were extracted and calibrated using the 
OPTEXT routines^ in conjunction with IRAF. The final 
flux calibrated spectra agreed well with the existing 5A spec- 
tra (Paper II) except in the case of BRI0952— 0115 where the 
signal-to-noise ratio was very poor and the relative flux dif- 
ferent by a factor of two. This spectrum has been excluded 
from further analysis. Though not of exceedingly high qual- 
ity the spectra were suitable for fitting the damped Lya 
candidates. 



2.2 Fitting the Damped Candidates 

The damped candidate system profiles were fit using the 
vpgtiQ Voigt profile fitting software. The programme re- 
quires three input spectra: the object, the errors, and the 
continuum. The OPTEXT routines used in the reduction 
created the first two and the continuum spectra were cre- 
ated by using continuum fits for the WHT spectra (Paper 
II) and extrapolating them to the echelle spectra as was 
done in Williger et al. 1994. The centroids of the Lya fea- 
tures were determined from narrow metal lines (e.g. OI, CII) 
and then the Lya lines were fit in a region 50- 100 A around 
the candidate feature. Generally, candidate damped systems 
seen in the low resolution spectra appeared much stronger 
than they actually are, due to blending of the dense Lya 
forest lines at low resolution. In all three QSOs the candi- 
date damped systems are seen at higher resolution to break 
up into multiple overlapping components. The complexity 
of the blend coupled with residual uncertainty in the exact 



placing of the continuum level causes some of the error esti- 
mates from VPFIT to be somewhat optimistic. In particular 
for BRI1108-0747 and BR1202-0725 numerical simulations 
suggest an error in log A^hi of a = 0.15 should be adopted. 
The fitted lines are summarised in table ^. 

a) BR1033-0327, (z em = 4.509, z absorptlon = 
4.15) The abs orption candidate a t z=4.15 had been previ- 
ously studied (Williger et al. 1994) with a 12 km s _1 echelle 
spectrum taken at CTIO. Their spectrum covered only the 



* These were written by Bob Carswell, Jack Baldwin, and Gerry 
Williger for the reduction of CTIO echelle data, 
t Written by R.F. Carswell, J.K. Webb, A.J. Cooke, and M.J. 
Irwin, see also Cooke (1994). 



blue wing of the system and from this the column density 
was estimated to be no greater than 3 x 10 20 atoms cm -1 . 
With spectral coverage of the entire absorber it is seen to be 
a complex system of at least 5 absorbers with a total column 
density of log N m = 20.15 ± 0.11 atoms cm -2 (1.4 x 10 20 ). 
In Paper II the column density was estimated to be log 
A?hi= 20.2. The redshifts of the Lya absorbers were deter- 
mined from 5 tentatively identified OI1302 lines in the CTIO 
spectrum. The signal-to-noise ratio in the order where the 
Ly/3 lines lie was too low to use in the fit. The spectrum 
with the best fit line profiles and ±l<r fits are shown as solid 
lines in figure 0(a). The error array and normalised con- 
tinuum are shown as dotted and dashed lines, respectively. 
The small spike at ~6270A is real. It also appears in the 5A 
resolution spectrum (Paper II). 

b) BRI1108-0747, {z em = 3.922, z Ma = 3.607) 
The absorber at z=3.61 is barely damped with a column 
density of log Nni= 20.33 ± 0.15 atoms cm -2 . In Paper II 
we estimated the column density to be log Nm= 20.2, so this 
system was not originally in the statistical sample. Several 
of the absorbers in the survey in Paper II have estimated 
column densities near the statistical sample threshold of log 
Nm= 20.3. We expect some to be confirmed above this value 
and some below. The profile fit and ±l<r fits are shown in 
figure |f](b). The redshift centroid was determined from a 
single strong CII 1334 line. 

c) BR1202-0725, (z em = 4.694, z dla = 4.383) 
The damped system in this QSO has a column density of 
log A m = 20.49 ± 0.15. It is the highest redshift damped 
Lya system known. The ESO spectrum is shown with the 
profile fit and zbla fits in figure |f|(c). It was also measured 



in a higher resolution spectrum taken at CTIO (Wampler 
et al. 1996), and includes over 10 components. Lu et al. 



(1996) have observed it at Keck with HIRES and measure 
a column density of log Nm = 20.6 ± 0.1. It was estimated 
in Paper II to have a column density of log Nm= 20.5. 
The status of the systems detected in the APM Damped 
Lya absorption survey is summarised in table ^| (It is an 
updated version of table 7 from paper II.) Those marked 
with an asterisk have column densities log Ahi> 20.3 atoms 
cm -2 and make up the statistical sample of high redshift 
absorbers used in the analysis. In Paper IV this data set is 
combined with previous surveys to study the evolution of the 
cosmological mass density of neutral gas at high redshift and 
the implications for galaxy formation theories are discussed. 
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Table 3. Status of Absorbers in APM Damped Lyman-a Absorption Survey 



QSO 




^TTLCIX 




z dla 


W rcst 
A 


log N m 
estimate actual 


BR 0019 - 1522 


2.97 


4.473 


4.528 


3.42 


7.6 


zu.u 












3.98 


12.3 


zu.o 












4.28 


8.0 






BRI0103 + 0032 


2.87 


4.383 


4.437 


4.23 


5.8 


ry.o 




BRI0151 - 0025 


2.74 


4.142 


4.194 










BRI0241 - 0146 


2.86 


4.002 


4.053 


3.41 


5.7 


ty.o 




BR 0245 - 0608 


2.96 


4.186 


4.238 










BR 0351 - 1034 


3.09 


4.297 


4.351 


3.62 


6.6 


1 o o 
ry.y 












4.14 


6.3 


1 o o 
ty.y 




BR 0401 - 1711 


2.82 


4.184 


4.236 










BR 0951 - 0450 


2.93 


4.315 


4.369 


3.84 


24.0 


*91 O 












4.20 


10.6 


*9fl Q 
ZU.O 




BRI0952 - 0115 


2.99 


4.372 


4.426 


4.01 


18.0 


ZU.O 




BRI1013 + 0035 


2.61 


4.351 


4.405 


3.10 


17.5 


*9n c 












3.73 


9.6 


90 9 

zu.z 












4.15 


7.4 


90 n 




BR 1033 - 0327 


2.91 


4.454 


4.509 


4.15 


9.6 


90 9 


on i ^ J. n 11 


BRI1050 - 0000 


2.83 


4.233 


4.286 










BRI1108 - 0747 


2.64 


3.873 


3.922 


2.79 


8.3 


90 1 












3.607 


9.0 


90 9 


ZV.oo ± U.lo 


BRI1110 + 0106 


2.58 


3.869 


3.918 


3.25 


5.2 


1 Q 7 

iy. / 












3.28 


6.2 


1 Q Q 

ry.y 




BRI1114 - 0822 


3.19 


4.440 


4.495 


3.91 


6.7 


1 Q Q 

ly.y 












4.25 


11.7 














4.45 


5.3 


1 Q 7 

iy. i 




BR 1144 - 0723 


Removed from 


sample. 










BR 1202 - 0725 


3.16 


4.637 


4.694 


3.20 


5.4 


1 Q 7 

iy. / 












3.38 


7.1 


90 n 
zu.u 












4.13 


7.8 


90 1 












4.383 


13.2 


90 K 


*9D /1Q 4- n K 


BRI1328 - 0433 


2.24 


4.165 


4.217 


3.08 


8.3 


20.1 




BRI1335 - 0417 


3.08 


4.342 


4.396 










BRI1346 - 0322 


2.65 


3.942 


3.992 


3.15 


6.8 


19.9 












3.36 


5.0 


19.7 












3.73 


10.0 


*20.3 




BRI1500 + 0824 


2.39 


3.894 


3.943 


2.80 


11.3 


*20.4 




GB 1508 + 5714 


2.73 


4.230 


4.283 










MG 1557 + 0313 


2.66 


3.842 


3.891 










GB 1745 + 6227 


2.47 


3.852 


3.901 










BR 2212 - 1626 


2.69 


3.940 


3.990 










BR 2237 - 0607 


2.96 


4.502 


4.558 


4.08 


11.5 


*20.4 




BR 2248 - 1242 


2.94 


4.109 


4.161 











* These absorbers are above the statistical sample threshold of TVhi > 2 x 10 20 atoms cm~ 2 . 
z m i„ = minimum redshift at which a DLA could be observed 
z em = emission redshift of the QSO 
z-max = 3000 km s^ 1 blueward of z em 

z dla = redshift at which a damped candidate was observed 
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3 THE HI COLUMN DENSITY 

DISTRIBUTION FOR LOG N m > 17.2 

3.1 Background 

The distribution of the HI column densities for QSO ab- 
sorption line systems has been investigated by several au- 
thors. Tytler (1987) found that the distribution may be rep- 
resented by a single power law f(N) oc N~ /3 over the range 
13.3 < log N m < 21.8 with = 1.51 ± 0.02. He argued that 
this was evidence for a single population. Fitting the higher 
and lower column density systems separately, Tytler found 
P = 1.61 ± 0.10 for log N m < 16.7 and (3 = 1.35 ± 0.07 for 
log Nyo> 17.2, very similar slopes particularly considering 
the inadequacies of the data set. Sargent et al. (1980) ar- 
gued that the forest lines and higher column density metal 
absorbers were distinct populations due to their different 
clustering properties. Sargent, Steidel & Boksenberg (1989) 
found a single power law fit the over the entire column den- 
sity range but noted that since little detailed data was avail- 
able for 15.5 < log JVhi< 17.2, the region where the HI be- 
comes optically thick, that this was not evidence for a single 
population. They also found j3 = 1.39 for log Ahi> 17.3. 
LWTLMH found the best fit for the high column density 
data (17.2 < log 7V H i< 21.8) gave p = 1.25. All of these 
values for /3 agree within the errors. Petitjean et al. (1993) 
reviewed these results and concluded that neither a single 
nor double power fit well for 13.3 < Ahi< 21.8 but that 
there clearly was a flattening of the distribution function 
around log JVhi~ 16. 

Looking at the damped Lya systems alone (20.3 < log 
N m < 21.8) LWTLMH found (3 = 1.73. Most of these investi- 
gators find no evidence for any substantial redshift evolution 
in the distribution function though LWT find an increased 
number of the highest column density systems at high red- 
shift. Into this quagmire we wade with the absorption sys- 
tems from the APM Damped Lya Survey [Paper II] and the 
results for the Lyman-limit system evolution from Paper I 
to attempt to better quantify what is happening with the 
high column density HI at high redshift. 



Table 4. Redshift and Absorption Distance Paths 



Data Set Az AX AX 

(q = 0) (go = 0.5) 



< z < 4.7 



APM Damped Lya Survey 


36.1 


162.3 


76.5 


WTSC + LWTLMH + LWT 


203.4 


602.1 


344.8 


Combined 


239.5 


764.4 


421.3 


2 > 3 








APM Damped Lya Survey 


30.5 


141.2 


65.5 


WTSC + LWTLMH + LWT 


13.0 


54.6 


26.7 


Combined 


43.5 


195.8 


92.2 



(Bahcall & Peebles 1969; cf. Tytler 1987). The value of go has 
little effect on the slope of the column density distribution 
function. In the analysis below we have utilised go = 0.5. 
To allow for redshift evolution, the distribution function is 
normally generalised as 

f{N,z) = kN- /3 (l + z)~<. (5) 

The damped Lya systems and candidates from the APM 
survey shown in table ^| have been combined with previ- 
ous surveys for damped Lya systems (WTSC, LWTLMH, 
LWT) which results in a data set of 366 QSOs yielding 44 
damped systems with log Ahi> 20.3 covering the redshift 
range 0.2 < z < 4.4. The total redshift path and absorp- 
tion distance covered by the surveys is shown in table ^. A 
maximum likelihood technique, described in Appendix A, 
has been employed to find values for and k to determine 
if a power law fit describes the HI column density distribu- 
tion for this sample. As already indicated a disadvantage of 
a power law model for the HI column density distribution 
of damped Lya absorbers, is the divergent nature of the 
integral mass contained in the systems. Since it is straight- 
forward to generalise the maximum likelihood method to 
alternative forms of the distribution, we explore in section 
3.6 an alternative parameterisation based on a gamma dis- 
tribution (cf. Pei & Fall 1995). 



3.2 A Power Law Distribution Function 

The form of the power law column density distribution func- 
tion f(N) typically used is 

-0 



f(N) = kN~ 



(2) 



f(N)dNdX is defined as the number of absorbers in an ab- 
sorption distance interval dX with HI column density N in 
the range N + dN . The absorption distance X is used to 
remove the redshift dependence in the sample and put ev- 
erything on a comoving coordinate scale since Az = 0.5 at 
redshift 2 is not the same at as a Az = 0.5 at redshift 4. If 
the population of absorbers is nonevolving (i.e. their num- 
ber density multiplied by their cross-section does not change 
with redshift), the absorption distance can be defined as 



X(z)= ={l + z){l + 2q a z)- 1/2 dz, 
Jo 

therefore 



X{z) 



|[(l + z)^ 2 -l 
![(l + *) 2 -l] 



if q = 0.5; 
if qo = 0. 



(3) 



(4) 



3.3 Results for a Single Power Law 

A single power law fit to the combined data set described 
above of the form in equation ^ with go = 0.5 results in 
= 1.74 ± 0.12 and logfc = 13.9 ± 2.3 with similar values 
for go = (P = 1.75, log k — 13.9). The quoted error for the 
normalisation constant k is large because for a small change 
in the value of P, k can change by 2 orders of magnitude. 
These results are in good agreement with the results found 
by LWTLMH (0 = 1.73 ± 0.29, log k= 13.63 ± 0.09) and 
are plotted for the entire data set in figure ^. The data is 
binned for display purposes only with the vertical error bars 
plotted at the mean column density for each bin. We will 
see in the next section that a single power law is not a good 
fit to the data. 



3.4 Cumulative HI Distribution 

As shown in Paper I for the Lyman-limit system evolution, 
the arbitrary binning of the data for presentation in differ- 
ential plots includes a subjective component that can mask 
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exactly what is happening in the underlying data. A cumu- 
lative distribution plot is far better at revealing the true 
nature of the distribution and this approach is examined 
now. The log of the cumulative number of damped Lya sys- 
tems detected versus log Nbi is plotted in figure 0(a). A 
point for the expected number of Lyman-limit systems that 
would be detected down to log Ahi= 17.2 is shown with a 
circled star. This is calculated by integrating the number 
density per unit redshift (N(z) = 0.27(1 + z) 1 ' 55 ) over the 
redshift path covered by the n QSOs in the DLA sample, 



LLS ea: p e cied = ^ / N(z)d Z = ^ J N {1 + z)"' dz 

= E / 0.27(1 + z) 1 - 55 dz. 



(6) 

It is obvious from figure 3(a) that a power law will not fit 
the entire column density range 17.2 < log JVhi< 22. A 
Kolmogorov-Smirnov (K-S) yields a probability of less than 
10~ 7 that the fit represents the underlying data set. In fig- 
ures ^(b-d) the same distribution is overplotted with single 
power law fits for different values of f3 that were fit to the 
graph by eye. (b) shows that (3 = 1.34 will fit from the Ly- 
man limit column density through the damped systems with 
log -/Vhi~ 21, a flatter slope than the canonical 1.5-1.7 range, 
(c) shows that f3 = 1.69 fits the damped distribution with 
20.3 < log -/V H i< 21.2 well but does not describe the high 
or low column density tails of the distribution, (d) shows 
a fit to the sharp drop off in numbers for damped systems 
with log Am Z 21.3. This can be expected from looking at 
the estimated column densities for the damped systems in 
table ^| or by looking at the spectra. There are not a lot 
of heavily damped systems. Clearly the results for a single 
power law fit depend critically on the range of column den- 
sities included. This characteristic can explain much of the 
variation in the results previously seen by various authors 
that were summarised in section 3.1. 



3.5 Redshift Evolution for 0.008 < z < 4.7 

To qualitatively study the redshift evolution of the column 
density distribution of the damped Lya systems the cumu- 
lative distribution shown in figure ^ has been split in half at 
redshift 2.5 and each set plotted individually (figure^). The 
damped Lya systems with z > 2.5 are shown by the solid 
line and the absorbers with z < 2.5 are shown by the dashed 
line. The higher redshift absorbers appear to have a slightly 
flatter slope up to log A r Hi=21 and then a sharper drop in 
the number of very high column density systems, though 
a K-S test shows that this is not a statistically significant 
difference. The evolution with redshift in the slope of the 
column density distribution is also apparent when looking 
at the differential f(N). LWT plotted this in 3 redshift bins 
z=[0.008,1.5], [1.5,2.5], and [2.5,3.5]. In the highest redshift 
bin there was a flattening of the column density distribution 
slope towards higher column densities. In figure [| we have 
plotted our combined data set with this same binning with 
the addition of one higher redshift bin z= [3.5,4.7]. The flat- 
tening of the distribution function towards higher column 



density systems in the z=[2.5,3.5] bin in the LWT data is no 
longer pronounced. The most striking feature is the steep- 
ness of the distribution in the highest redshift bin. It is not 
just steeper due to a decrease in the highest column den- 
sity systems (log A r Hi>21), but there is also an increase in 
the number of lower column density systems relative to the 
other bins. Even if 15-20% of the candidate systems with 
log A r Hi~20.3 turn out not to be damped when observed at 
higher resolution, as we expect, this result still holds. 



3.6 Results for a Gamma-Distribution 

There are two strong motivating factors to find an alterna- 
tive model for describing the HI column density distribution. 
First, as shown in section 3.4, there is direct evidence for 
an apparent variation in the power law slope as a function 
of iVm- This implies that a higher order functional form 
other than a power law is needed to describe the column 
density distribution. Second, as noted in the introduction, 
with a power law model the integral mass contained within 
damped Lya systems is divergent for realistic values of (5. 
This in turn means that it is impossible to assign a formal 
upper limit to any estimate of the neutral gas content of the 
early Universe. Consequently, following Pei and Fall (1995), 
we have chosen to model the data with a gamma distribution 
of the form 



f(N, 



(f«/N*)(N/N t 



-P -N/N, 



(7) 



where /» is the characteristic number of absorbing systems 
at the column density N t , and N* is a parameter defining 
the turnover, or 'knee', in the number distribution. Both /* 
and N* may in general vary with redshift but for the moment 
we treat them as constants. This functional form is similar 
to the Schechter luminosity function (Schechter 1976). For 
N « N* the gamma function tends to the same form as 
the single power law, f(N) cx A r_ ' 3 ; whilst for N > N*, the 
exponential term begins to dominate. 

We can understand how a gamma function might pro- 
vide a better description of the damped Lya data by con- 
sidering the differential logarithmic slope which is given by 



(8) 



dlogf(N,z) 
dlogN 

As the column density approaches N, the slope begins to 
steepen and rapidly turns over at higher column densities, 
qualitatively similar to what we observe in figures ^| and ^| 
Furthermore, the integral HI over the column density dis- 
tribution (cf. equation [j]) for N m i n « N, is now given by 
/*JV»r(2 — 7) where r denotes the standard gamma function. 
This function is bounded if 7 < 2. 

The maximum likelihood technique outlined in Ap- 
pendix A can readily be modified to incorporate this form. 
We note that the likelihood solution can be found over a 
two-dimensional grid of pairs of values of JV* and f3, since 
the constant /* can be directly computed using the con- 
straint 



f(N,z)dzdN 



(9) 



where m is the total number of observed systems. This is 
computationally much less intensive than doing a 3-D grid 
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search. The results of a single functional fit to the entire 
dataset are log N, = 21.63 ± 0.35, f3 = 1.48 ± 0.30, and 
/* = 1.77 x 10~ 2 . The log-likelihood function results with 
confidence contours are shown in figure ^(a). The best fit is 
overplotted on the differential form of f(N) in figure |^(b) 
and on the cumulative distribution in figure Be). (The single 
power law form of f(N) was shown fitted to the same data in 
figures and ^|). The differential form of the plots (figures [| 
and |^(b)) show little difference between the single power 
law and T-distribution fits. When displayed with the cumu- 
lative number of absorbers in figure |6|(c), the T-distribution 
now clearly fits the entire data set with column densities log 
JVhi> 20.3. If the expected number of Lyman- limit systems 
are included in the fit, the results are log N* = 21.36 ±0.15, 
P = 1.16 ± 0.15, and /» = 4.43 x 10~ 2 . This also provides a 
reasonable fit to the data as shown in figure (?](a-c). 



4 NUMBER DENSITY EVOLUTION WITH 
REDSHIFT 

Differential evolution in the number density of damped Lya 
absorbers has been described by LWT and Wolfe et al. 
(1995). While the change in number density per unit redshift 
is consistent with no intrinsic evolution of the absorbers over 
the range < z < 3.5, they find that the systems with log 
./Vhi> 21 disappear at a much faster rate from z=3.5 to z=0 
than does the population of damped absorbers as a whole. 
We now examine the redshift evolution of the damped Lya 
absorbers in our combined data set by determining the num- 
ber density of absorbers per unit redshift, dN/dz = N(z). 
In a standard Friedmann Universe for absorbers with cross 
section tvRo 2 and number density <E>o per unit comoving vol- 
ume 



N{z) = $oitRo 2 cHq\1 + z)(l + 2q Q z)- 1/2 . 



(10) 



It is customary to represent the number density as a power 
law of the form 



N(z) = N (1 + Z y 



(11) 



where iVo = QqhR^cHq . This yields 7 = 1 for qo — and 7 
= 1/2 for qo — 1/2 for the case of no evolution with redshift 
in the product of the number den sity and cross section of 
the absorbers ( gargent et al. 198C ). 

A maximum likelihood fit to the data yields N(z) = 
0.04(1 + z) 1 °' B which is consistent with no intrinsic evo- 
lution even though the value of 7 is similar to that found 
for the Lyman limit systems where evolution is detected 
at a significant level (Paper I; Stengler-Larrea et al. 1995). 
The log-likelihood function for 7 and No with >68.3% and 
>95.5% confidence contours is plotted in figure []. We also 
find redshift evolution in the higher column density systems 
but with a decline in N(z) for z>3.5. These results are dis- 
played in figure |^. The entire data set is plotted as dashed 
lines with the above fit. The results for only the absorbers 
with log N(HI)> 21 are shown as solid lines. Figure |ic| shows 
HI column density versus redshift, and the paucity of ab- 
sorbers with log Nm> 21 at z > 4 is apparent. 



5 CONCLUSIONS 

Three QSOs from APM survey have been observed at 0.8A 
resolution. Two have damped systems with confirmed HI 
column densities of Nm> 10 20 ' 3 atoms cm~ 2 , with a third 
absorber falling just below this threshold. We have dis- 
covered the highest redshift damped Lya absorber known 
at z=4.383 in QSO BR1202-0725. The two systems with 
N m > 10 20 ' 3 atoms cm , and remaining nine candidate 
damped absorbers from the APM survey have been com- 
bined with data from previous surveys to study the column 
density distribution and number density evolution for ab- 
sorbers with iVm > 17.2. If the HI column density distri- 
bution function is fit with a power law, f(N) = kN -13 , 
we find evidence for breaks in the power law, flattening 
for 17.2 < log iVm ~ 21 and steepening for log Nm ~ 
21.2. The column density distribution function for the data 
with logiVHi>20.3 is better fit with the T-distribution form 
/(AO = {f t /N t ){N/N t )~ l3 exp{-N/N t ) with log N, = 
21.63 ± 0.35, P = 1.48 ± 0.30, and /» = 1.77 x 10~ 2 . 

For the number density evolution of the damped ab- 
sorbers (log Nm > 20.3) over the redshift range 0.008 < z 
< 4.7 we find the best fit of a single power law form for 
N(z) = N (l + z)~< yields 7 = 1.3 ± 0.5 and N = .04t° 3 . 
This is consistent with no intrinsic evolution in the absorbers 
even though the value of 7 is similar to that found for the 
Lyman limit systems where evolution is detected at a signifi- 
cant level. Evolution is evident in the highest column density 
absorbers with the incidence of systems with log N(HI)>21 
decreasing for z > 3.5. 
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APPENDIX A: MAXIMUM LIKELIHOOD 
ANALYSIS 

Using equation |^ for the column density distribution func- 
tion, the damped Lya absorbers will be found randomly 
distributed according to this function along the QSO line- 
of-sight in N — z space. If the space is divided into m cells 
each of volume Sv, the expected number of points in cell i 
is given by 

4n = f(N,z) t 6v. (Al) 
The probability of observing Xi points in cell i is 

p(x,) = e-^J. (A2) 

The likelihood function for QSOj taking the product over 
all the cells is then 



(A3) 



If the volume of each cell Sv becomes very small such that 
there is either 1 or points in each cell, 



I 1 ' 
10, 



if DLA detected; 
if none detected, 



then the likelihood can be rewritten separating out the terms 
for full and empty cells. For m = g empty cells + p full cells 



IK* IT 



IK* IT 



(A4) 



i — 1 j — 1 i — 1 j — 1 

Taking the log of the likelihood function gives 

m p 

lo S L 3 = ^2 ~<t>i + ^2 ln 4>3 
i=l 3=1 
m p 

= 5^-/(JV,*)ifo + j^In/(J\T,z)j+pln&> 



3 = 1 



(cf. Schechter & Press 1976). Ignoring the constant terms, 
in the limit where Sv — *• this becomes 



log L/, 



N m , 
N min 

N min 



f(N,z)dNdz + J2\nf(N,z) j 
j=i 

p 

kN~ {l + zY'dNdz + ^ln[fcAT _/3 (l + z) 



3=1 

(A6) 

To get the overall log lik elihood for n QSOs we evaluate the 
integrals in equation A6 and additively combine the log L's 
resulting in 



log 



l = E[ (1 -^r +7) (( i+ ^ 1+7 -( i+ ^ 1+7 ) + 



pin k 



i=l 
Pi 



+ g(-/31nAr 3 .+ 7 ln(l + 4j) 



3=1 



(A7) 



where pi is the number of detected DLAs in QSO; and Nmin 
is the minimum column density. 



(A5) 
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Figure 1. The profile fit and ±1<t fits to the damped Lya absorbers are shown as solid lines. The error arrays are shown as dotted lines. 
The components are listed in table 2[ (a) The HI absorption complex in BR1033— 0327 at z=4.15 is a system of at least 5 absorbers with 
a total column density of log Nm= 20.15 ± 0.11 atoms cm -2 , (b) The damped Lya absorber in BRI1108— 0747 at z=3.607 with log 
iV HI = 20.33 ± 0.15 atoms cm -2 , (c) The damped Lya absorber in BR1202-0725 at z=4.383. This is the highest redshift damped Lya 
absorber known. The central damped component is shown with log 7Vjji=20.49±0.15. 



Figure 2. A single power law form of the column density distribution function, f(N) = kN @, fit to the entire damped Lya sample 
from the APM Damped Lya Survey, WTSC, LWTLMH, and LWT. The parameters of the fit are /3 = 1.74 and logfc = 13.9. 



Figure 3. (a) The cumulative distribution for 17.2 < log Nm < 22. The stepped line is the data for all the damped Lya systems in the 
data set. The circled point is the number of Lyman limit systems that would be expected given the redshift path covered in the damped 
Lya surveys. In (b)-(d) the same distribution is overplotted with single power law fits for different values of /3 that were fit to the graph 
by eye. (b) shows that /3 = 1.34 will fit from the Lyman limit column density through the damped systems with log TVjji ~ 21, a flatter 
slope than the canonical /3 1.5 — 1.7 range, (c) shows that j3 = 1.69 fits the systems with 20.3 < log TVhi < 21.2 well but does not 
describe the high or low column density tails of the distribution, (d) shows a fit to the sharp drop off in numbers of damped systems 
with log TVjji Z 21. This is evident from just looking at the estimated column densities for the damped systems in table jj| or by looking 
at the spectra. There are not a lot of heavily damped systems. 
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Figure 4. The cumulative HI column density distribution with the sample split in half at z = 2.5. The damped Lya systems with 
z > 2.5 are shown by the solid line and the absorbers with z < 2.5 are shown by the dashed line. The higher redshift absorbers appear 
to have a slightly flatter slope up to log TVhi = 21 and then a sharper drop in the number of very high column density systems, though a 
K-S test shows that this is not a statistically significant difference. 



Figure 5. The log column density distribution function /(TV) vs. the log column density TVhi plotted over 4 redshift ranges, z=[0.008, 
1.5], [1.5, 2.5], [2.5, 3.5], and [3.5, 4.7], all binned in the column density ranges log 7V H i=[20.3, 20.5], [20.5, 21.0], and [21.0, 21.8]. The 
gradual flattening of the distribution function from redshift z=0 to z=3.5 is evident. The most striking feature is the steepness of the 
distribution in the highest redshift bin. It is not just steeper due to a decrease in the highest column density systems (log TVjji>21), but 
there is also an increase in the number of lower column density systems. 



Figure 6. (a)The log-likelihood function for the T-distribution form of the column density distribution function, /(TV, z) = 
(f*/N*)(N/N*)~Pexp(-N/N*). The >68.3%, >95.5%, and >99.7% confidence contours are plotted for log TV, and (3. The best fit 
values are log TV* = 21.63, /3 = 1.48, and /» = 1.77 X 10~ 2 which is solved for analytically, (b) The T-distribution function, 
/(TV,z) = (/* /N*)(N/N*)~Pexp(—N/N*), overplotted on the differential form of the column density distribution. The fit is to the 
entire data set including the surveys from Paper II, WTSC, LWTLMH, and LWT. The parameters for the fit are log TV, = 21.63, 
fi = 1.48, and /* = 1.77 X 10 -2 . (c) The cumulative distribution for 17.2 < log TVjji < 22 as shown in figure ^(a) is now overplotted 
with the T-distribution form of the column density distribution function. The dashed line now clearly fits the entire distribution for log 
A r jjj>20.3. The circled point is again the number of Lyman limit systems that would be expected given the redshift path covered in the 
damped Lya surveys. This is not included in the fit. 



Figure 7. This figure shows the same data plotted in figure |fj but the fit includes the expected number of Lyman limit systems, 
given the redshift path surveyed, (a) The log-likelihood function for the T-distribution form of the column density distribution function, 
/(TV,z) = {f*/N*)(N/N*)-Pexp(-N/N*), The >68.3%, >95.5%, and >99.7% confidence contours are plotted for log TV* and /3. The 
best fit values are log TV, = 21.36, f) = 1.16, and /» = 4.43 X 10 — 2 which is solved for analytically, (b) The T-distribution function, 
/(TV, z) = (/* /N*)(N/N*)~@exp(—N/N t ), overplotted on the differential form of the column density distribution. The fit is to the entire 
data set including the surveys from Paper II, WTSC, LWTLMH, and LWT. The parameters for the fit are log TV, = 21.36, /3 = 1.16, 
and /» = 4.43 X 10 -2 . (c) The cumulative distribution for 17.2 < log TVjji < 22 overplotted with the T-distribution form of the column 
density distribution function. The circled point is again the number of Lyman limit systems that would be expected given the redshift 
path covered in the damped Lya surveys. 



Figure 8. The >68.3% and >95.5% confidence contours for the log-likelihood function are plotted for the number density per unit 
redshift of the damped absorbers. The best fit for a single power law form TV(z) = TVq(1 + z) 7 yields 7 = 1.3 ± 0.5 and TVo = .04 j[ Q2 over 
the redshift range 0.008 < z < 4.7. 



Figure 9. The number density of DLA per unit redshift, TV(z), vs. z (absorption). The dashed bins show N(z) for all the damped systems 
and the solid bins for systems with N(HI)> 10 21 atoms cm~ 2 . A single power law fit of N(z)= .04(1 + z) 1,3 is overplotted. This is 
consistent with no intrinsic evolution in the absorbers even though the value of 7 is similar to that found for the Lyman limit systems 
where evolution is detected at a significant level (Paper I). 



Figure 10. The HI column density of the damped Lya absorbers is plotted versus absorption redshift. The paucity of absorbers with 
log TVhi> 21 at z > 4 is apparent. 
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